Exact joint density-current probability function for the asymmetric exclusion process 
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We study the asymmetric exclusion process with open boundaries and derive the exact form of 
the joint probability function for the occupation number and the current through the system. We 
further consider the thermodynamic limit, showing that the resulting distribution is non-Gaussian 
and that the density fluctuations have a discontinuity at the continuous phase transition, while the 
current fluctuations are continuous. The derivations are performed by using the standard operator 
algebraic approach, and by the introduction of new operators satisfying a modified version of the 
original algebra. 
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Due to the lack of a general theory of non-equilibrium 
steady-states, a lot of the research in this area focuses 
on the study of simple models. Of special interest are 
stochastic interacting particle models Q, Q, and one 
hopes that many of the interesting qualitative character- 
istics of these simplified models are generic for a larger 
class of systems. A system that has received a lot of 
the attention is the partially asymmetric exclusion pro- 
cess (PASEP). Firstly this model is non-trivial, display- 
ing steady-state phase transitions, yet simple enough to 
be integrable 0, fl EL0 , and further it maps onto cer- 
tain growth models |J it models traffic flow 8], and 
it is believed to describe the large scale dynamics of 
the noisy burgers equation |jj 0, [n| and the KPZ 
equation 0, There has been much progress in 

the analytical treatment of the PASEP, giving rise to 
host of exact results describing it's steady state prop- 
erties H H H H [3 H [H At the heart of the 
distinction between non-equilibrium and equilibrium sys- 
tems lies the ability of non-equilibrium systems to carry 
currents. So far, the results concerning the currents in 
different special cases of the PASEP [H Q d 03 are 
mainly for systems with periodic boundaries or infinite 
geometries with special initial conditions. In this letter 
we consider a finite system with open boundaries, but 
specialize the treatment to the one dimensional asym- 
metric exclusion process (ASEP). It consists of a lattice 
of size L, with site label I running from left to right. 
Every site on the lattice can be occupied by no more 
than one particle. Given that the right neighboring site 
of an occupied site is empty, the occupying particle will 
jump to the empty site with a rate 1. If the first site 
on the lattice is unoccupied, particles are injected at this 
boundary with rate a. Further given that we have a par- 
ticle at the last site of the lattice, it is ejected with the 
probability rate (3. No further transitions are allowed. 
We will here limit our considerations to the case where 
we can view the boundary rates as deriving from parti- 
cle reservoirs. We therefore take < a — pi c f t < 1 and 



deviations, open systems, stationary non-equilibrium 



< /3 = 1 - pnght < 1, where pi cft and p righ t are the 
particle densities of the reservoirs. This model has been 
exactly solved (see Figure ^for the phase diagram) in 
the sense that the steady-state probability of any given 
microscopic configuration can (in principle) be calculated 
by applying a given set of algebraic rules. Even so, these 
calculations quickly become very cumbersome as the sys- 
tem size is increased. Thus we here wish to extract gen- 
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FIG. 1: Phase diagram of the one dimensional exclusion pro- 
cess. The dashed line indicates the first order transition line, 
while the dash-dotted lines indicates the continuous transition 
lines. 

eral information about the system directly from the alge- 
braic rules, without explicitly calculating the microscopic 
weights. Since the algebraic rules are instrumental to 
our later development we here give a very brief recap on 
their definition. The starting point is to represent any 
microscopic configurations in terms of a string of non- 
commuting operators D and E, corresponding to a par- 
ticle and a hole respectively. It can then be shown that 
the steady-state probability distribution can be written 
in terms of this operator string and two vectors, (a\ and 
|/3), according to 



P ss ({ ni }) = {Zf )- l {a\X{m)X{n 2 ) 



■X(n L )\p). (I) 



Here the operator X{ni) equals D if there is a particle 
at site I (ni — 1), and E if site / is unoccupied (n; = 0). 



The state independent factor Z" = (a\(D + E) L \(3) en 
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sures the proper normalization. For (JTJ to hold true, the 
operators and vectors must further satisfy the algebraic 
rules 

C^DE = E + D, (a\E = -(a\, D\p) = i|/3), (2) 

a [3 

where we have implicitly assumed that the normaliza- 
tions of the vectors \a) and (/3| are such that (a\/3) = 1. 
The algebraic rules (0) are now all that is needed to cal- 
culate P^({ni}), resulting in a polynomial of degree L 
in 1/a and 1/(3. 

In moving from a microscopic to a macroscopic view 
of the system we will here concentrate on the average 
density and current throughout the bulk. We derive the 
exact joint probability function for the average bulk cur- 
rent and density for any system size. First we define the 
total activity within the system as the number of bulk 
bonds that can facilitate a transition of a particle in the 
immediate future, i.e. the total effective bulk transition 
rate. The bulk current is then defined as the activity di- 
vided by the system size. For any given state the activity 
equals the number of pairs of neighboring sites that has a 
particle to the left and a hole to the right. To get a han- 
dle on the activity, J, of a microscopic configuration of 
N particles we choose to represent such a configuration 
by a sequence of J objects of the form D p 'E h i , possibly 
padded with E's to the left and .D's to the right. Doing 
this we can write any microscopic steady state measure 
as 

P ss {{ni}) = 

(Z^)- l {a\E ho (D Pl E hl ) ■ ■ ■ (D p 'E hj )D Po \p), 

by appropriately choosing the numbers {pj , hj } and J. 
It further follows that the above expression is unique if 
ho,Po > and the rest satisfy hj,pj > 1. We can now in 
principle calculate the joint probability distribution for 
N and J by summing the above over all hj 's and pj 's con- 
sistent with a specific number of particles (Y)j—Q Pj = N) 

and a given system size (N + J2j=o hj = L). Choosing 
to enforce these constraints with contour-integral repre- 
sentations of the Kronecker delta, the expression for the 
joint particle- activity probability function can be written 



dzdz 
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P a P(N J) = — 

^ Z L J (27TZ) 2 Z N+1 E L-N+1 



C z ,Cz 



N 



L-N 



(z — 0){z — a) 



J2^DrJ2(zE) h ) |/3). (3) 



0=1 



h=l 



Here C z (C s ) is a directed contour that encircle the pole 
at the origin of the complex z (z) plane once in the pos- 
itive direction, with \z\ < (3 (\z\ < a). The first step 
toward explicitly calculating Q is through considering 



the properties of the operators ^Z p {zD) p and ^ h {zE) h . 
Surprisingly it turns out [lij that a slight modification 
of the above operators 



D' ={l-(z + z))D N Y;\zD) h , 

h=0 

a r L-N-l 

E> d ^(l-(z + z))E £ (zE) p , 

p=0 



satisfy the "relaxed" operator algebra 

D'E' = D' + E' + 0{z N ,z L - N ). 

The relaxed eigenvectors and eigenvalues are further 
given by 

D'\P) = 0(z N ), (a\E> = ±j{a\ + 0(z L - N ), 

p' a' 

with 



, def a — z dof (3 - z 

a = ; -r, [3 = 



l-{z + z)' 



1- z + z) 



The fact that these eigenvalues are complex is of no con- 
cern since we consider only finite polynomials in the in- 
verse eigenvalues. Any result is thus uniquely extendable 
into the complex plane through analytic continuation. 
We can now rewrite Q in terms of the primed oper- 
ators, and start using the relaxed operator algebra to 
transform the expression. The result of any such manip- 
ulation would be the same, up to terms of order z N and 
z L ~ N , as if the operator algebra would have been exact. 
The extra terms have no effect under the contour inte- 
gral in (PJ) since the poles at the origins are both of order 
equal or lower than N and L — N. (The case for J = is 
trivial.) Thus, using the relaxed algebra to perform any 
manipulation within is equivalent to using an exact 
algebra. Therefore we can write 



Pf(N,J) 



a/3 



dzd 



z 



{2m) 2 z 



2 ^N+l-J~L-N+l-J 



7a'&' 



(z - (3){z - or)(l - (z + z)) 



2J ' 



(4) 



This expression is the main result of this letter, and 
since all quantities in it are known exactly, it yields 
both the exact finite system size form of P^(N, J), as 
well as the asymptotic form in the large system size 
limit. By first calculating the generating functional 
G(M) = T,^=oV L zf El (which has recently been de- 
rived in a different manner in .20]), the above further 
yields a double contour inte gral expression for the gener- 
ating functional of N and J |Tjj • Below we present exact 
and asymptotic results for P^f (N, J) . 
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Finite systems: The integral in J2J is easily calculated 
with the help Cauchy's integral theorem. All we need to 
do is to calculate the coefficient of the term proportional 
to (zz)" 1 in the Laurent-series expansion of the integrand 
in (2) . For J = we have Z 



a'/3' 



1 and thus 



Pf(N,0) 



(l/(3f(l/a) L - 



N 



which is obviously correct since the inactive state must 
have L—N empty sites followed by N filled sites. For J > 
1 we use the exact form of the normalization function 



Z^ = > ' A 



L 

E- 

i=i 



i-LJ 



1 1 

E a kM- 
fe=0 h 



A L ,i = 



l(2L-l-l)\ 
L\{L-l)\ 



to write 



^(^) = f§E^E E E 

L j=l k=0 c=0 d=Q 



j L-N-JN-J 



Gh,c(ot)Gj-k,d(P)H2J-j,L-N-J-c,N-J-di (5) 



with 



Get 
k,c 



1 



vc+fc+l 



H 



K,d,e = 



c J a 
K-l + d + e\/d + e 
d+e )\ e 



Through the above we now have the exact form of the 
sought-after joint probability function for any system 
size. The form is illustrated in Figure |2 

The thermodynamics limit: We here return to 
Using the asymptotic form of the normalizing function 
given in 0, we perform a steepest descent calculation 
to get the asymptotic results. We consider the different 
phases individually. Due to the particle-hole symmetry 



Pf{N,J) 



a (L — N, J) it is only necessary to ex- 
plicitly consider the case a < (3. 

First turning to the maximal-current phase we con- 
sider and drop all pre-factors that are independent of 
N and J (this will be done throughout) to write 



Pf(N,J) 



dzdz 



1 



J 3 / 2 / (2m) 2 Z N-J+1 E L-N-J+1 

1 1 

(1 - {z + z))^- 1 (2a - 1 + z - z) 2 {2[3 - 1 - {z - z)) 2 ' 

The asymptotic behavior of these integrals is in principle 
straight forward to calculate. In practice though, it turns 
out to be quite cumbersome since one has to determine 
which of the saddle points and lower order poles give the 
dominant contributions. We can shortcut this through 
only considering the asymptotic form in some finite re- 
gion around the peak of the distribution. From we 




a = .5, B = .5 
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FIG. 2: Each row contains a surface and a contour plot of 
the exact probability distribution for the values of a and (5 
indicated, and with p = N/L and j = J/L. The first three 
rows illustrate the behavior of the probability distribution as 
the system goes along the line of a = /3 through the critical 
point at a = /3 = .5, while the last three graphs illustrate the 
behavior as the system goes through the first order transition 
at a = f3 — 0.25. Overlaid in the contour plots (dashed line) 
is the curve j — p(l — p) which defines the set of possible 
asymptotic average values of p and j throughout the systems 
different phases (not at the first order transition line). The 
system size is L = 40. 



know that the average density and current is a and (3 
independent. Thus, the lower order poles cannot dictate 
the asymptotic behavior around the peak value of the 
probability distribution, and instead this must be set by 
the saddle points 

z*=p-j, z* = l-p-j, p = N/L, j = J/L. 

A saddle-point approximation thus results in 



Pf(pJ) 



1 



f3( p -j)P-3(i- p -jy-P-i 



(G) 



where we for simplicity have dropped all the sub- 
dominant pre-factors. Even though the extent of the re- 
gion of validity of 10 is unknown, it should be pointed 
out that the size of this region is finite (as long as the 
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system is away from any phase boundaries) and indepen- 
dent of system size. In the first row of Figure |3 we show 
the resulting dominating asymptotic plots. 




FIG. 3: Fhe two rows display a surface and a contour plot 
of the leading behavior of the asymptotic probability distri- 
bution. The calculations were performed at the injection and 
ejection rates indicated and at a system of size L = 40 (to 
make the result comparable to Figure 

Now turning to the low-current phases we have 



be a finite jump in the connected density-density correla- 
tor. Comparing this to equilibrium systems, these tran- 
sitions correspond to proper second order transitions. 

In conclusion we note that it would be very interest- 
ing to examine if the same "trick" of introducing a re- 
laxed algebra could somehow be applied to the PASEP. 
This especially since this model interpolates between a 
equilibrium and non-equilibrium steady state. It would 
further be interesting to derive the full asymptotic form 
of the probability distribution since sufficiently close to 
a phase transition, any finite system will reach a point 
at which the region of validity of the above asymptotic 
forms shrink to the size of the typical fluctuations. In 
this region the system crosses over to a situation where 
the asymptotic fluctuations are governed by the tails ex- 
cluded in the above development. 
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from the Merton College Oxford Domus scholarship fund, 
and the Royal Swedish Academy of Science. 



dzd 



z 



1 



P L (P,j) ' i Z N-J+1 Z L-N-J+1 

C z ,Cz 

2a - 1 + z - z 1 



0-a-(z-z)(a- z) J+1 (l - z ~ a) J+1 ' 

The same arguments as applied in the maximal-current 
phase now gives us the asymptotic probability distribu- 
tion around the peak. Again it is the saddle points 

p - 3 s -* 1 - p - i 

z = (1 — a), z = a 

P 1 - P 

that dominate. The resulting dominant form is 



Pf(N,J) 



P p {l-pf- p 



-p(1 - a)P jV(p - j)P-i{l - p - jf-P-i 



(7) 



The above result, valid for the low-density phase, is di- 
rectly transferable to the high-density phase through the 
use of the particle hole symmetry mentioned above. A 
realization of the asymptotically dominating part in the 
low-density phase is shown in the second row of Figure^ 
It is clear from the asymptotic forms that the prob- 
ability distribution is non-Gaussian in all phases. This 
is consistent with the view that long-range correlations 
are a generic feature of non-equilibrium systems |2l| . In 
general, as a phase-transition line is approached the bor- 
der of the region of validity of the asymptotic forms © 
and must approach the peak. It should be further 
pointed out that as the continuum transition is passed, 
the asymptotic forms JfJJl and 0) indicates that there will 
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